

%% parameters
load SMM
fprintf('verif poids=1');
p0 =[ p01_SMM;p02_SMM];
mu_c = [mu_c1_SMM; mu_c2_SMM];
%sig_c = sig_c_SMM;
sig_eps_a = [sig_eps_a1_SMM; sig_eps_a2_SMM];


    [f1, avfracup1,    med1,    interq1, kur1, f2, avfracup2,    med2,    interq2, kur2]=geNCalvoPlus_SMM(p0, mu_c, sig_eps_a, weight_j);
    
    
    
    